Collapse of Small-Scale Density Perturbations during Preheating in Single Field 

Inflation 
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After cosmic inflation and before the transition to radiation domination, the cosmic energy density 
may have been dominated during an extended period by an oscillating massive scalar condensate. 
We show that during this period, sub-Hubble scale perturbations are subject to a metric preheating 
instability in the narrow resonance regime. We analyze in detail both, quadratic and quartic po- 
tentials. The instability leads to the growth of density perturbations which in many cases become 
non-linear already before the beginning of a radiation dominated Universe. This is particularly the 
case when requiring a phenomenologically preferred low reheat temperature. These early structures 
may lead to the emission of gravitational waves and the production of primordial black holes. Fur- 
thermore, it is not clear if they could modify the prediction of linear curvature perturbations on 
very large scales. 



I. INTRODUCTION 



The Universe having passed through an early inflation- 
ary epoch is by now a widely accepted paradigm for an 
explanation of apparent cosmic large-scale homogeneity 
as well as the perturbations observed in the cosmic mi- 
crowave background radiation (CMBR). One of the earli- 
est incarnations of inflation was given in Ref. [l[ where a 
massive scalar field with potential V{4>) = m 2 (/> 2 /2, start- 
ing from super-Planck values, slowly rolls down towards 
its minimum thereby driving a period of accelerated ex- 
pansion. Though many alternate scenarios of inflation 
have been discussed over the years, the above scenario of 
"chaotic inflation" remains one of the most successful in 
terms of reproducing the observed perturbations [2|. In 
particular, the predicted spectral index n s ~ 0.96 in this 
model is consistent with that inferred from the WMAP- 
7 analysis of the CMBR inhomogeneities, which gives 
n s ~ 0.963 ± 0.012 (68% CL) 0. 

After inflation has ended, cosmic energy density is 
dominated by an almost homogeneous oscillating infla- 
ton field. A transition towards a radiation dominated 
Universe, as required to exist at the latest shortly be- 
fore the epoch of Big Bang nucleosynthesis (BBN) which 
starts at cosmic temperature T ~ 1 McV, is usually as- 
sumed to occur by the perturbative decay of massive in- 
flaton particles into radiation. This epoch is referred to 
as the epoch of reheating. Parametric resonance phe- 
nomena, inducing energy transfer to other scalars cou- 
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pled to the inflaton, referred to as preheating, though 
rich in their physics, have been shown to fail to provide a 
complete transition towards a radiation dominated Uni- 
verse (see, e.g. 0, Q for a review of preheating mod- 
els). In case the inflaton would decay immediately af- 
ter the end of chaotic inflation, and assuming m 2 (j) 2 /2 
inflation, the cosmic reheat temperature would be ex- 
ceedingly large, i.e. T r h ~ 10 16 GeV. Such a large T r h 
is phenomenologically problematic, as likely leading to 
the re-generation of unwanted relics, such as magnetic 
monopoles at a GUT transition, or massive gravitinos in 
supergravity 0, Modern constraints from BBN on 
unstable gravitinos constrain this reheating temperature 
down to T Th < 10 7 GeV [HljJ. It is thus possible (and 
desirable) that between the end of chaotic inflation and 
reheating cosmic energy density is dominated for a very 
long expansion period by an almost homogeneous infla- 
ton condensate. For example, for T r h as low as ~ 1 MeV 
the Universe would have expanded an immense factor 
~ 10 25 between the end of inflation and reheating. It is 
such a cosmic era which is the topic of the current paper. 

An inflaton oscillatory phase after inflation has re- 
ceived some attention with respect to the evolution of 
super-Hubble size perturbations [IiTfl3l ]. It has been 
shown that, though one finds parametric resonance in the 
m 2 cj> 2 /2 case, super-Hubble size perturbations thought to 
be responsible for the anisotropies in the CMBR stay con- 
stant in amplitude, as is the case in a radiation dominated 
Universe. It was therefore concluded that in single-field 
models "metric preheating" does not lead to non-trivial 
changes in the evolution of perturbations. However, as 
we will see in this paper (see also Ref. [Hj]), parametric 
resonance occurs even for sub-Hubble size regions, i.e. for 
all modes with physical wave vectors fc p h ys < VSHm (for 
a quadratic potential), where H is the Hubble constant. 
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In fact, these modes lie in the narrow resonance regime 
and the red-shifting is such that any mode which had 
once been resonant will always remain resonant. Due to 
this resonance the curvature perturbations of sub- Hubble 
modes stay constant rather than decreasing. One can 
readily show that this implies that the density contrast 
grows with scale factor Spk/p oc a implying non-linear 
cosmic structure formation after an expansion factor of 
~ 10 6 . This fact has so far received little attention in the 
literature. 

Whether the Universe indeed entered an era with non- 
linear structures, already before the epoch of reheat- 
ing, depends on the reheat temperature and the scale 
of inflation. Attempting to keep the reheat tempera- 
ture low enough to avoid overproduction of gravitinos, 
i.e. T rh < 3 x 10 7 - 3 x 10 8 GeV @, [HI implies that es- 
sentially all large-field inflationary models are followed 
by a phase of non-linear structure formation. Figure Q] 
provides an illustration of this fact assuming a quadratic 
potential during the oscillating phase. It is seen that 
even for H cn a as low as 10 7 GeV non-linearities typically 
develop unless T r h > 10 8 GeV. Therefore, a nonlinear cos- 
mic phase before reheating appears to be a generic pre- 
diction. Note, however, that these considerations assume 
that after inflation energy density is dominated by a sin- 
gle oscillating scalar field. 

These findings may not be surprising since it is known 
that the effective equation of state during such a phase 
is that of a Universe filled with pressure-less dust. One 
would thus expect the inflaton condensate to be gravita- 
tionally unstable, which indeed, has been shown in the 
Newtonian limit [l4], EH and in an averaged relativistic 
limit 0,113 • However, our analysis and the effect that 
we report go beyond this Newtonian analogue, because 
this instability is found to be intrinsically some form of 
metric preheating. And, indeed, we find that a similar 
effect also takes place in the case of a quartic inflaton 
potential, for which the effective equation of state dur- 
ing oscillations is that of radiation. In this latter case, 
the resonance is (and remains) sharply peaked around 
the physical scale associated to the Hubble scale at the 
end of inflation, but the growth of the density contrast 
is exponential. 

The outline of the paper is as follows. In Section |TT] 
the metric preheating of scalar perturbations is analyzed 
in detail. The primordial power spectrum is calculated 
numerically over a wide range of scales. In Section IIII1 
we examine in detail the growth of density perturbations 
in a quadratic inflaton potential, while in Section IIV1 we 
discuss the case of a quartic inflaton potential. Conclu- 
sions are drawn in Section W\ 



II. GENERAL CONSIDERATIONS 

The evolution of scalar (density) perturbations is con- 
trolled by the so-called Mukhanov variable. If matter is 
described by a scalar field (as it is the case during in- 
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FIG. 1: Contour plot of the r.m.s. density fluctuation Spk/p 
on the scale corresponding to the horizon at the end of infla- 
tion as measured at the onset of the radiation era in the T r h 
(reheating teniperature)-i/ en d (Hubble constant at the end of 
inflation). The solid line that is labelled "1" delineates pa- 
rameter space where non-linear structures form between the 
end of inflation and reheating (the region below this line) 
from parameter space where the Universe remains in the lin- 
ear regime of density perturbations, i.e. Sp^/p < 1, before 
reheating (the blue region above the line). The normalization 
for the curvature perturbations on the Hubble scale at the 
end of inflation has been assumed to be Vc, ~ 10 -11 in accor- 
dance with Fig. [21 The two dotted lines indicate where the 
power spectrum of the fluctuations is larger and smaller than 
one or, equivalently, how the solid line would move when the 
normalization of the curvature perturbations on the Hubble 
scale after inflation is higher, Vc, = 10™ 9 (line labeled 0.01), or 
lower, V( = 10~ 13 (line labeled 100). The solid red line shows 
approximate upper limits on T r h in order to avoid gravitino 
overproduction. Note that the figure is based on Eq. (|25|) . and 
makes the fairly generic assumptions that energy density after 
inflation is dominated by a single oscillating scalar field with 
potential m 2 </> 2 /2, and that <j> < M pl at the end of inflation. 



flation and preheating), then the equation of motion of 
the Fourier transform of Mukhanov's variable is given 
by Pill 



k 2 - 



v k = . 



(1) 



In this article a prime denotes a derivative with respect to 
conformal time while a dot means derivative with respect 
to cosmic time. The quantity k is the comoving wave 
number and t\ = —H/H 2 is the first slow-roll parameter. 
The quantity H = a/a is the Hubble parameter and a 
denotes the Friedman-Lemaitre-Robertson- Walker scale 
factor. Mukhanov's variable is related to the curvature 
perturbation, 



K Vk 



2 a^/ei 



(2) 
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where k = 87r/m P1 , mp\ being the (non-reduced) Planck 
mass. The importance of the curvature perturbation lies 
in the fact that, on large scales, it is a conserved quan- 
tity [iH and, therefore, it can be used to propagate the 
inflationary spectrum from the end of inflation to the 
post-inflationary era. As a consequence of the above re- 
lation, the spectrum of Ck can be expressed as 



for instance due to particular symmetries of the underly- 
ing theory. Different choices of coefficients in the above 
equation ([6]) correspond to different models. With the 
three terms displayed in Eq. ©, one can describe large 
field models, small field models and hybrid inflation (at 
least in the inflationary valley). In the next section, we 
consider the massive case. 
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Vk 
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In order to calculate V((k), one needs to integrate 
Eq. (TTJ), which requires the knowledge of the initial con- 
ditions for the mode function vu- Since, at the begin- 
ning of inflation, all the modes of astrophysical interest 
today were much smaller than the Hubble radius, the ini- 
tial condition are chosen to be the Bunch-Davis vacuum 
which amounts to 



III. THE QUADRATIC CASE 
A. Density Perturbations 

In this section, we investigate a model with a simple 
quadratic potential, 



lim Vk 

fc/(aiJ)-H- oo 



-ikr) 



(4) 



1 



where r\ denotes conformal time. 

The above formulation is well-adapted to the calcula- 
tion of the power spectrum during inflation. However, in 
order to study the metric perturbations during the phase 
where the inflaton field oscillates it is more convenient 
to work in terms of cosmic time and to use the rescaled 
variable defined by Vk = a}/ 2 Vk- Then, Eq. (TTJ) takes the 
form [13! 



Vk + 



„2 



d 2 V 




+ 2k 



<f> dV 
~H~d4 



v k = 0. 



(5) 



This formulation is particularly convenient as all the 
terms in the above equation remain non-singular even 
when the field oscillates at the bottom of its poten- 
tial, that is to say when <f>, <f> and e% periodically van- 
ish [Ull. 

To proceed further, one needs to choose the potential 
V(<p). The form of the inflaton potential is presently 
unknown (although it is constrained, see Refs. (23 - [26| ). 
A convenient way to treat the most general case is to 
consider the following potential 



V{4>) = V a ± imV + 



(6) 



where the dots refer to higher order terms. As it is well- 
known, this expansion can be difficult to control, most 
notably so for large field models where, during inflation, 
the vacuum expectation value of the field is large com- 
pared to the Planck mass. However, one can always as- 
sume that some terms in the above expansion are absent, 



V{(j)) = -m A 4> 



2i2 



(7) 



We have in mind a large field inflationary model, but as 
discussed above, this quadratic potential has a more gen- 
eral applicability when it comes to discussing the physics 
of inflaton oscillations after the actual era of inflation. 
It may also describe the physics of curvaton oscillations. 
We will therefore keep wherever needed, the dependence 
of our results on the mass m (and Hubble scale H en d 
at the end of inflation) , although we will use as fiducial 
values those corresponding to a large field inflationary 
model. 

In large field models, the mass m is fixed by the nor- 
malization of the amplitude of density fluctuations to 
that measured by COBE and WMAP. Standard calcu- 
lations lead to 



rms-PS 

m ~ mpi 

2iV* + 1 T 



(8) 



where iV» is the number of e-folds before the end of in- 
flation at which the scale of astrophysical interest left 
the Hubble radius during inflation (we take iV* = 50) 
and where Qrms-Ps/T - 6 x 10 -6 . This gives m ~ 
1.4 x l(T 6 mpi. 

The power spectrum 7\(fc) of density fluctuations at 
the end of inflation is calculated numerically in the fol- 
lowing, via the integration of Eq. ([1} throughout the in- 
flationary era. Although one may find excellent approxi- 
mations to this power spectrum on CMBR scales through 
the use of standard slow-roll calculations, we indeed find 
that up to second order in the slow-roll parameters, these 
methods cannot predict reliably the power spectrum on 
the shorter scales that we are interested in. More specif- 
ically, one finds to first order in the slow-roll parame- 
ter (23 
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1 - 2 (C + 1) ei - Ce 2 + (-2ei - e 2 ) In 
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P(W = 5- 

7reim P1 

while, at second order, it is given by [13] 

T c (k) = H \ (l-2(C + l)ei-Ce 2 + I 2C' Z + 2C-, , 
7reim P1 [ " V 2 



(9) 



5 ej + C - C 



7^ 
12 



7 eie 2 



C 2 ^ 2 



He 2 



2 



e 2 e 3 + [-2ei - e 2 + 2 (2C + 1) e 2 + [2C - 1) e l£2 + Ce 2 - C*e 2 e 3 ] 



x In 



(4e 2 + 2eie 2 + el - e 2 e 3 ) In 2 



(10) 



In these expressions, e 2 and e 3 are slow-roll parame- 
ters belonging to the Hubble flow hierarchy [27l l28j 
e n+ i = dlne n /(diV) (N being the number of e-folds) 
and C is a numerical constant C ~ —0.7296. The quan- 
tity fc* is the so-called pivot scale the location of which 
is arbitrary but is usually chosen in the middle of the 
range of scales probed by the CMBR. These approxima- 
tions fail to reproduce the power spectrum on small scales 
because the lever arm is then too important for the accu- 
racy of the slow-roll calculation. This is rather unusual, 
as the slow-roll method generally provides an excellent 
approximation. This break-down of the slow-roll calcula- 
tion on small scales is clearly seen in Fig. [2] which com- 
pares the full numerical calculation of the power spec- 
trum with the above approximations for a large field 
model quadratic inflationary model. As one might ex- 
pect, second-order slow-roll calculations are more accu- 
rate on small scales than the first order (which, among 
others justifies the techniques that have been developed 
for this purpose [29|, [ 



Evolution of Density Perturbations after Inflation 



We are now interested in the behavior of perturbations 
after the end of inflation. On time scales t 3> H~^ d (the 
subscript "end" denotes hereafter quantities evaluated at 
the end of inflation), the inflaton vacuum expectation 
value evolves approximately as 
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FIG. 2: Primordial power spectrum of curvature perturba- 
tions (solid line) for the large field model V(4>) = m 2 (j) 2 /2. 
Here iV e xit (fc) — -/Vend is In [ 

Across (fe)/oend], with across (k) the 
scale factor at Hubble radius crossing of mode k during in- 
flation (i.e. k/a CIOSS = 2nH), and a en d the scale factor at 
the end of inflation. The dotted and dashed lines respec- 
tively represents the first and second order slow-roll results 
as given by Eqs. © and JTDJ with e\ ~ l/[2(iV» + 1/2)], 
£2 = £3 = l/(A r » + 1/2) and N, ~ 50. The slow-roll result 
is accurate in the CMBR region but breaks down on smaller 
scales where the predicted power spectra become negative. 
The upturn is related to modes that have not crossed out the 
Hubble radius and have remained in their vacuum state, see 
text for details. 



(Qend \ 3/2 
-^-J sin (mi) , (11) 



for a quadratic inflaton potential. Let us recall that 
cn d ~ topi/(2y / 7t) and fiend — mj \[% in the particular 
case of a large field quadratic inflationary model. Plug- 
ging Eq. (TTTT) into Eq. ((5$ , it is easy to see that the dom- 
inant term of the last four gravitational terms is the last 
one as it scales as 1/a 3 / 2 while the others decay as 1/a 3 . 
Neglecting the other three terms one finds an equation 



of motion of the form 
d 2 v k 



dz 2 



fc2 fTTM f a cnd\ 3 / 2 ,„ v' 

-3-9 - V 6re<rWd I cos(2z) 



1 + '., .-, - V'i/'-'Vn,) ( — — ) cos(2z) Vk=0, 

where we have defined z = mt + 7r/4. This equation is 
similar to a Mathieu equation 



d 2 v k 
dz 2 



+ [A k - 2gcos(2z)]v fe = 



(13) 
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FIG. 3: Hubble radius £ H (red solid line), lower bound of 
resonance band £ c (red dashed line) and physical wavelengths 
(blue and green dotted lines), in Planck units, versus the e- 
fold number during the final stage of inflation and the infiaton 
oscillation era for a large field m 2 2 /2 model. The reheating 
temperature is T r h = f0 7 GeV hence the number of e-folds 
during the infiaton oscillation era is JV r j, — 25.6. The inset 
shows the detailed behaviors of £ H and l Q at the transition 
between inflation and infiaton oscillation. 



with 



the resonance band during inflation (blue lines in Fig. [3]) 
and those that have entered it during the infiaton oscilla- 
tion era (green lines in Fig. [3]). Of course, the preheating 
instability is only operative after inflation and before re- 
heating. 

Having identified the scales affected by the resonance, 
we now discuss their amplification. For the first instabil- 
ity band, the Floquet index is given by fx = q/2. Then 
the mode function evolves as Vk oc e MZ . However, in 
the present case, we have a time-dependent Floquet in- 
dex, for which the corresponding solution can be written 
as fil 



Vk oc exp 



fidz 



,3/2 



(17) 



Note that, taking the integral of the time-dependent /i 
does not correspond to the WKB approximation, but 
rather is an educated guess confirmed by numerical simu- 
lation, see below. This implies that Vk oc a and therefore 
(k is constant according to its definition Eq. Note 
also that the curvature perturbation usually remains con- 
stant only on super-Hubble scales while decaying on sub- 
Hubble scales, unless there is gravitational instability. 

We are now interested in the behavior of the fractional 
mass-energy density perturbation, (5^ (77) = Spk/p, where 
p is the background energy density of the scalar field. The 
perturbed Einstein equation implies that 



A k = 1 



/ 6k , /a en d\ 3 / 2 

TT~ <Pend 



(14) 
(15) 



Since q <C 1, we are in the narrow resonance regime. The 
first instability band is given by 1 — q < Ak < 1 + q which 
amounts to 



< - < VSHm . 
a 



(16) 



From the above inequality, we note the emergence of a 
new characteristic spatial length given by £ c = 1/VaS, 
which sets the lower bound to the range of resonant spa- 
tial scales. During inflation, in a large field model, £ c is 
almost a constant and £ c <~ £ H , with £ H = H^ 1 . After 
the end of inflation, during the infiaton oscillation era, 
with H cx a~ 3 / 2 one has £ c cx a 3 / 4 and, therefore, £ c 
grows more slowly than the Hubble radius £ H cx a 3 / 2 . 
The situation is summarized in Fig. [3] where both scales 
are displayed (solid and dashed red lines). Their evo- 
lution has been computed by an exact numerical inte- 
gration of the equations of motion in the case of a large 
field model. Since the physical wavelength of a comoving 
Fourier mode, A oc a, grows faster than £ c , the scales may 
cross £ c during inflation or during the matter-dominated 
era and thereby enter the instability band. Furthermore, 
a mode that has entered this band will always remain in 
it, at least until reheating is completed. One can thus 
distinguish two types of modes: those that have entered 



Sh = — - 



2 H 2 



H ' 



(18) 



where is the Bardeen potential jig, bL0] 
potential is related to the quantity Cfe by 



Cfe 



3 l+p/p 



k ■ 



The Bardeen 



(19) 



where p denotes the pressure. For constant p/p, if 
is constant so is Qk- This is a generic result on large 
scales in one fluid models. If $k evolves as a power-law 
in the scale factor (also a generic result), 3>fc oc a s , then 
i/ _1< f>fc = s$fc and consequently, Eq. (fT9"j) indicates that 
Cfc follows the same power-law behavior. Therefore, the 
constancy of in the resonance band implies that 
also remains constant. This is confirmed by a numerical 
integration of the equations of motion, see Fig. |4j Then 
Cfe = 5$fe/3, as the pressure vanishes on average during 
the preheating stage. Consequently, one may express 5k 
as a function of Oc as follows 



2 ( k 2 

Oh = 

5 



H 2 



3 Cfe- 



(20) 



This relation remains valid in both the large scale and 
the small scale limit (with respect to the Hubble radius) 
provided the modes are in the instability band (i.e. Cfe is 
constant). Given that aH oc a -1 / 2 in a infiaton oscilla- 
tion dominated Universe, one then finds that 8k grows as 
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FIG. 4: Evolution of the density perturbations Spk/p, and 
(frfc towards the end of inflation and during the ensuing period 
of metric preheating after inflation. The notation is as in the 
text with fcphys = k/a; N denotes the e-fold number of the 
scale factor. The reduced Planck mass Mpi = mpi/\/87r is 
used in this plot. 

a(t) on sub-Hubble scales but remains constant on super- 
Hubble scales. This behavior is also confirmed by the 
numerical integration whose results are shown in Fig. 3) 
The preheating instability for all sub-Hubble modes with 
k/a < \/3Hm may thus be understood as the gravita- 
tional instability of the oscillating scalar condensate. 

Let us summarize our main findings. We have found 
that all modes in the instability band and outside the 
Hubble radius are such that and Sk are constant. 
This is of course the case for modes on scales probed by 
the CMBR and, therefore, one recovers the fact that the 
inflationary power spectrum propagates through the re- 
heating stage without being affected. However, for those 
modes which are in the instability band and inside the 
Hubble radius, one finds that is constant but that Sk 
grows as a(t). Therefore, and contrary to the standard 
lore, the preheating stage can affect the behavior of the 
perturbations even if there is only one non-interacting 
fluid. 



2. Spectrum of Density Perturbations at Reheating 

From the above considerations, we see that the ampli- 
fication described above occurs only for a range of scales, 
k € [k-min, fcmax] in terms of comoving wavenumber. We 
now discuss this point in more detail and specify fc m ; n and 
fcmax- Henceforth, we introduce the short-hand notation: 

* = ^4— ' (2i) 

where the label "end" refers as before to the end of 
inflation. The smallest unstable spatial scale is that 



which enters the resonance band at the time of comple- 
tion of reheating. It corresponds to the green dotted 
line labelled "fc max " in Fig. [3j Interestingly enough, this 
scale has never been super-Hubble. From the condition 
fcmax/flrh = \/3iJ r hTO, where Hrh is the Hubble parame- 
ter at the end of the reheating epoch or at the beginning 
of the radiation dominated phase, one obtains 

/ m \ 1/2 

fe max = 1.37 x 10 3 

\ 1.4 x 10~ b mpi ) 

X (lO^GevJ (jO^GevJ ' (22) 

The last equation assumes <?* = 230 degrees of freedom 
at reheating. All modes with k < fc max have entered the 
resonance band during inflation or during the preheating 
era. 

We have also seen that a necessary condition for the 
growth of Sk on a given scale is that the corresponding 
mode lies within the Hubble radius at some time before 
the reheating epoch. As a consequence, the largest spa- 
tial scale that is unstable with respect to the growth of 
6k is that which re-enters the Hubble radius at the end 
of the pre- reheating epoch, i.e. k min ~ a r h-ff r h- This cor- 
responds to the blue dotted line labeled "fc m j n " in Fig. [3] 
and 

(rp \ 2/3 / TT \ - 1 / 3 

iPSv) (hSSv) (23 » 

As the density contrast increases in direct proportion to 
the scale factor for modes in the range k £ [fc m j n , fc max ], 
non-linear structures may eventually form. The ampli- 
tude of the density contrast at the end of the preheating 
epoch is given by the density contrast at the end of infla- 
tion and by the e-fold number of growth while the mode 
is sub-horizon sized and in the resonance band. At this 
point, one should distinguish the modes that enter the 
resonance band during inflation (the dotted blue lines in 
Fig. [3]) and those that have entered the band during the 
preheating era but that have never exited the Hubble 
radius (the dotted green lines in Fig. [3]). 

We first address the former ones, for which k 6 
[&min, fcend]) with fc cn d = 1. Given that, in the resonance 
band, the Fourier density contrast remains constant on 
super-Hubble scales, the behavior of 5k on sub-Hubble 
scale can be approximated as 

5k = 6l c — = ^ d — , (24) 

where the label "he" refers to the epoch of Hubble radius 
crossing of mode k and <5™ d is the density contrast calcu- 
lated at the end of inflation, see Fig. [5] Let us precise the 
value of ah c : in order to match both asymptotic behav- 
iors of Sk in the super-horizon, namely Sk — > 6£fc/5 and in 
the sub-horizon regimes, namely Sk — > [k/(aH)} 2 2^/5, 
one defines Ohc as that at which k — \f3ahcHhc- This 
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implies that a nc = 3a cn dfc 2 ■ Now, using Eq. (|20|) . one 
can write 6% nd = (6/5)^/2^^)7^. The n, the power 
spectrum of 8k at the time of reheating is evaluated as 



T)i-h/i \ a I rend I z I "rh 

27T V «hc 



2.82 x 10 lo fc 4 



H, 



end 



10 13 GeV 



io- n 

4/3 



-8/3 



10 7 GeV / 



(25) 



The variance of the density fluctuation enclosing the scale 
defined by fc is directly expressed at reheating as 



<^ 2 



din A:' Vl h (k') 



(26) 



and the condition (S 2 )k — 1 thus roughly amounts to 



1. This singles out a specific spatial scale (cor- 
responding to wavenumber fc n i), such that scales with 
k < k n \ have not had sufficient time to become non- 
linear, while those with k > k n \ have become non-linear 
by reheating. Straightforward manipulations lead to 



4.50 x 10" 



An 



10 13 GcV 



10" 11 

-1/3 



-1/4 



\ 10 7 GeV/ 



2/3 



(27) 



Let us now address the unstable modes that have al- 
ways been sub- Hubble, namely k G [fc nd, fcmax]- For 
these modes, Mukhanov's variable is still given by the 
vacuum solution, Eq. (j4)), until it enters the instability 
band. This leads to 



-ifc(r/-77 end ) 



(^f ) cr 1 , (28) 



where e\ ~ 3/2 in the preheating era. In order to derive 
this expression, we have set ei = 1 at the end of infla- 
tion and we have assumed that Mukhanov's variable has 
remained constant through the (instantaneous) transi- 
tion from inflation to the preheating era; this is expected 
insofar as the modes are well sub-Hubble at that time, 
k ^> aH, hence Mukhanov's variable oscillates freely. 
The quantity £™ d represents the value of at the end 
of inflation or at the onset of the oscillatory phase. Us- 
ing Eq. ([TO]) with p — 0, one finds £fc = since (k and 
must share the same power-law behavior oc a -1 (as 
discussed above), hence i? _1 $fc = Then, Eq. (fl"8|) 

leads to 



8k = —- 



2 fc 2 



3a 2 H 2 



-ik(r)-r} end ) 



~Sfe 



end 



C(fc)e 



-ikrf 



(29) 



The density contrast then oscillates with a constant am- 
plitude until the mode enters the resonance band, in 
which 8k — C(fc)a/a rc where a rc denotes the scale fac- 
tor at the entry into the resonance band. This latter 
is given by a rc = k 4 / (9m 2 al nd H 2 nd ) which implies that 
flrh/flrc = (kmax/k) ■ The power spectrum of density 
fluctuations at the time of preheating can now be com- 
puted as 



v?{k) = 



Vl nA {k) 



(30) 



The quantity 'p™ d is the power spectrum at the end of 
inflation for the modes that have never exited the Hubble 
radius, 



■pen, 



k 2 f -ffend 
( k ) = — 

tt \ mpi 



(31) 



We can check that it gives the correct order of magnitude 
by comparing with the exact result in Fig.[5]for the case of 
a large field model. Inserting this expression in Eq. (|30p , 
we see that the spectrum now decreases as k 2 while it 
was proportional to k in the other range of wavenumbers 
[fcmin > fcend] that we discussed previously. 

Finally, for those modes that have never entered the 
resonance band and never crossed out the Hubble ra- 
dius, that is to say k > /c max , the density contrast always 
behaves according to Eq. (j2"9")l and, therefore, the corre- 
sponding power spectrum is given by 



vT{k) 



27 



(32) 



This spectrum is proportional to fc 6 . Of course, these 
modes are still in a vacuum state and therefore, despite 
the fact that the vacuum fluctuations are large (which is 
likely to be linked to the pathological behavior of the the- 
ory in the ultra-violet regime) , there is no energy density 
associated to this phenomenon. One way to see this is 
to compute the evolution of the modes, and in particular 
the Bogoliubov coefficient of the negative frequency 
WKB branch at late times, on sub-horizon scales (as- 
suming the vacuum corresponds to the positive frequency 
branch). For these vacuum modes, j3k = because the 
WKB approximation has always remained valid. The 
energy density stored in these modes, which is given by 



d8pk ^ fc 4 
dlnfc ~ 27r 2 a' 



-Wk\ 



(33) 



consequently vanish. 

The situation is summarized in Fig. [2] where the spec- 
trum Vs(k) is displayed in the different regimes studied 
before. The striking feature of this plot is of course that, 
for fc > fc n i, the power spectrum is non- linear at the time 
of reheating; let us stress that of course, the plotted val- 
ues only represent the extrapolation of the linear calcu- 
lation. 
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k/KndHend) 



FIG. 5: Power Spectrum 7 , J h (fc), computed at reheating as 
given by Eqs. (|25[) . (|29|l and (|32l) . In the range k £E [A<min, ^end] ) 
the spectrum is evaluated numerically, see Fig. [2] The small 
discontinuity observed at fc en d is just a numerical artefact, the 
actual power spectrum being of course continuous. 



The non-linearity of modes with k > k cn a is rather 
intriguing, as these modes have never crossed out the 
Hubble radius, hence one might question their classical 
nature. However, it is rather straightforward to show 
that for these modes, \Pk\ ~ 1, hence energy density is in- 
deed stored on those scales. From this point of view, the 
non-linearities observed in the power spectrum should a 
priori be taken seriously and viewed as a sign that struc- 
ture formation has started on those scales. In order to 
assess the nature of these modes, one should apply the 
discussion of the so-called quantum to classical transi- 
tion of modes [2l], [3^] to the present case. However, this 
question is beyond the scope of this paper. 



3. Associated Mass Scales 

Finally, it is interesting to quantify the range of masses 
associated with the collapsing structures. In terms of the 
wavenumber of the perturbation, the enclosed mass can 
be written as 



Af=|(2.) 3 ^fc- 3 . 



H 3 

n end 



(34) 



According to the previous considerations, the mass corre- 
sponding to maximum growth of the density fluctuations 
corresponds to the scale fc en( j. Using the above expres- 
sion, straightforward manipulations lead to the following 
expression 



Me„d - 3.3 x 10 J 



H, 



end 



f ;ev) 



(35) 



The minimal mass, associated to k max is given by 
M( fcmax ) * 1.3x10- (^) 

-3/2 



(36) 



1.4 x 10- 6 mpi, 
The maximum mass is associated to the scale k n \, i.e. 

-2 

g- (37) 



M n] ^ :-!.(, ]()'"( JjL-V 7 (- 7>h 



Vl0 7 GeV 



Quite interestingly, the abundance of primordial black 
holes with mass < 10 9 g is constrained by early Uni- 
verse cosmology. Such black holes indeed evaporate be- 
fore BBN has started and the black hole temperature is 
high enough to allow the production of copious amounts 
of lightest supcrsymmetric particles, in excess of the clo- 
sure density [33j , as well as gravitinos and other late time 
decaying massive particles, whose decay products would 
ruin the success of BBN 34, 35] . If such black holes leave 
Planck-mass relics at their final point of evaporation, ad- 
ditional constraints can be obtained from closure density 
arguments 36]. For black hole masses > 10 9 g, stringent 
constraints can be derived from the impact of the evapo- 
rated radiation on BBN, see [37} for a detailed discussion 
of this and other observational constraints. 



B. Tensor modes 

In the previous section, we have studied the amplifica- 
tion of density perturbations during the reheating epoch. 
Naturally, one may wonder whether tensor modes can be 
amplified through a similar instability. Their Fourier am- 
plitude hk — Hk /a obeys the equation [38| 



/'-fc 



" 'fc 2 -- )m* = o. 

a 



(38) 



The effective potential is thus much simpler than in the 
case of density perturbations as it only involves the scale 
factor and its time derivatives up to second order. Defin- 
ing p,k = a 1 / 2 /ifc and working in terms of the cosmic time, 
one obtains 
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H 2 Ufc - . 



(39) 



Using the explicit form of the inflaton field given by 
Eq. (fTT|) . the equation of motion takes the form 



dz 2 



k 2 



3k o fa. 



end 



COS (2z) 



where z = mt. Therefore, one also obtains a Mathieu 
equation with 

k 2 



A k = 



3k 2 
_ 16 d 



(41) 
(42) 
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These equations should be compared to Eqs. (|T4|) 
and (fl"5|) . We are still in the regime q -C 1 but the cru- 
cial point is that there is no term 1 in the expression of 
Ak in Eq. (|4"Tj) . This means that, for a given scale, the 
condition 1 — q < < 1 + q cannot be maintained since 
Ak does not scale as q does. Therefore the amplification 
does not occur for gravitational waves. 



IV. THE QUARTIC CASE 

A. A quartic potential V(4>) = A</!> 4 /4 

In the previous section, we have studied a large field 
model with a quadratic potential. This corresponds to 
a pretty generic situation since, close to its minimum, 
the inflaton potential can always be approximated by a 
parabola. However, it might also be that, for some un- 
specified reason (generically, the presence of a particu- 
lar symmetry), the quadratic term in the Taylor expan- 
sion (J6j) is absent. Then, one has to repeat the previous 
analysis in the case where the potential is given by 



(43) 



In the following, we find that growth also occurs for this 
case. This result is particularly interesting, as it shows 
that the growth of density perturbations goes beyond the 
simple analogue of Newtonian collapse in a dust Universe. 
Indeed, in the case of A</) 4 /4, the effective equation of 
state is that of radiation. 

As before, we assume that the potential is approxi- 
mated by the above quartic form in the preheating era 
without making any specific assumption on the actual po- 
tential during the inflationary era. We note that quartic 
inflation can be normalized to COBE anisotropies with 
A ~ 10~ 12 , but that at the same time such an inflation- 
ary potential is disfavored by detailed analysis of CMB 
data [13. 

For a quartic potential, the evolution of the inflaton 
field around its minimum can be expressed in the large 
time limit t > H m \ as 0, i, M S3 

4>{rf) = $6 e nd (~~) Cn ( X ~ ^end! A=j , (44) 

where x = a en d ^end's/A??, f] being the conformal time and 
cn denotes the Jacobi cosine function [4l], [42[ . This solu- 
tion should be inserted in Eq. ([5]) . Then, keeping only the 
dominant terms, the equation that controls the evolution 
of Mukhanov's variable reduces to 



a 2 



— + 3A</> 2 nd 



•^end 



'V2 



v k = 0. 



(45) 



It is more convenient to change back to the variable Vk 
and to express this equation in terms of the variable x 
defined above: 



Vk 



dx 2 



k 2 



*ond a cnd^ 



3cn^ 



I 

71 



Vk = . 



(46) 

As emphasized in Refs. 0, [H, Ho|, this equation is the 
so-called Lame equation for which the Floquet chart 
is exactly known. There is amplification if 3/2 < 
k 2 /(<PL d a 2 cnd \) < V3. Using H 2 nd ~ the 
above inequality can be cast into the following form 



mpi 



2v^< 



< k < 



'end 



'3V3" 



1/2 



mpi 



(47) 



For the particular case of quartic large field infla- 
tion, </>end/"tpi — l/V''") hence the lower and upper 
bound of the above band reduce respectively to 3/2 and 
(3V3/2) 1 /2. 

For the above instability, the maximum growth occurs 
with a Floquet index \i ~ 0.036 at fc 2 /(</>g nd a 2 nd A) ~ 
1.615 (i.e. fc_~ 1.56 for a large field quartic inflation- 
ary model) |5|. The Floquet index is constant, therefore 
growth becomes exponential, unlike what has been found 
for the quadratic case. 

Several remarks are in order at this point. Firstly, the 
resonance band corresponds to a fixed narrow range of co- 
moving scales, contrary to the quadratic case for which 
the band widens with time. Secondly, the modes that are 
subject to parametric resonance are of spatial scale com- 
parable with the Hubble radius at the end of inflation. 
Thirdly, the effect found here is similar to that investi- 
gated in Ref. [J, Q except that it concerns Mukhanov's 
variable rather than the perturbed inflaton field or the 
Fourier amplitude of a scalar field x coupled to the infla- 
ton as <7 2 (/> 2 % 2 /2, where g is the coupling constant. The 
case of density perturbations can be viewed as a special 
case with a specific value of g. But, of course, the under- 
lying physics is different in the sense that metric fluctu- 
ations are different from an external scalar field coupled 

to (j). 

In order to check the effect discussed above, we have 
numerically integrated the equation of motion for a large 
field inflationary quartic potential, see Fig. [5] The nu- 
merical analysis confirms the previous analytical analysis, 
although we find that the instability band is shifted by 
12% to higher values in comparison with Eq. (14T1) ; this 
effect is likely due to our assumption of instantaneous 
transition from inflation to a (radiation dominated) pre- 
heating era in the analytical calculations. In detail, one 
finds the band at 1.684 < k < 1.807. When the mode is 
outside the resonance band (left panel in Fig. [BJ), the per- 
turbations remain constant while when it is inside (right 
panel in Fig. [6]), we observe a violent amplification of 
Cfe and 5 k - 

To summarize, the quartic case presents substantial 
differences in comparison with the quadratic case. There 
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FIG. 6: Evolution of the density perturbations Spk/p, (k and towards the end of quartic inflation and during the ensuing 
period of metric preheating after inflation. The dimensionless parameter is A ~ 10 -12 . The red solid line represents the Hubble 
radius while the blue dotted line represents the physical wavelength of the Fourier mode. Left panel: the Fourier mode is 
characterized by k ~ 1.427 and lies outside the resonance band. Right panel: k ~ 1.738 and the mode is inside the instability 
band as the exponential growth indicates. 



exists a narrow instability band at scales of order of the 
Hubble radius at the end of inflation, of fixed width (in 
terms of comoving wavenumber) and with associated ex- 
ponential growth. 

B. The mixed quadratic-quartic case 

The purely quartic case arises only if the quadratic 
term is absent in the potential. Here we briefly address 
the more general case, in which both quadratic and quar- 
tic terms are present: 

V{cf>) = \m 2 tf + . (48) 

At large vacuum expectation values of the inflaton, the 
quartic term always dominates while at small values of </>, 
the quadratic term is the most important one. For which 
value of 4> these terms are of equal importance depends 
on the two free parameters to and A. If the potential 
were purely quartic, then inflation would stops at <f> cn d — 
mpi/y/ft- For this value, the two terms are equal if m — 
\J A/(27r)mpi. It is therefore convenient to parameterize 
the potential (|48|) by A and the dimensionless parameter 
T defined by 




The regime where T ^> 1 reduces to the quadratic case 
while the case T <C 1 is equivalent to the quartic case. 
The most interesting situation is of course when T is of 
order one. 

In order to study this case, we have numerically inte- 
grated the equations of motion, assuming that the poten- 
tial form remains the same during inflation and the pre- 
heating era. As a first example, we have taken A = 10 -12 , 



that is to say the same value as in the previous section 
for numerical applications and T = 0.05. We have also 
chosen k = 1.738, i.e. a value which is inside the reso- 
nance band in the pure quartic case, see right panel in 
Fig. |H The result is displayed in Fig. [7] (left panel). Sev- 
eral comments are in order. Firstly, one notices that the 
Hubble radius and the characteristic scale l c no longer 
coincide. This is of course due to the fact that the mass 
to is no longer COBE normalized. In fact, in this model, 
the end of inflation is delayed because at <fi = mpi/y/w, 
the quadratic term still plays a role. Regarding the per- 
turbations themselves, despite the fact that the Fourier 
mode is inside the quartic resonance band of Eq. (|47p . 
there is no amplification. It is of course due to the fact 
that, for small vacuum expectation values of the field, 
the dominant term is the quadratic one. We also notice 
that 5k does not grow as a either. This is because, since 
the end of inflation has been postponed and the value of 
£ c increased (in comparison with the corresponding value 
in the pure quadratic case), the time at which the mode 
crosses l c is also delayed (and therefore lies outside of 
the plot). 

It is clear that if we decrease the value of T, and keep 
the same value for k, we should observe the growth de- 
scribed in the last section in the pure quartic case. This 
is indeed the case as demonstrated in Fig. [7] (right panel) 
where T = 0.003. Just after the end of inflation, the 
Fourier mode is inside the pure quartic resonance band 
and the perturbations grow. But as the amplitude of the 
inflaton field decreases, the influence of the quadratic 
term becomes more and more important. As a con- 
sequence, the exponential growth stops at some point. 
Then, the amplitude of the fluctuations remain constant 
until the mode exceeds the characteristic scale £ c asso- 
ciated to the quadratic part of the potential. As before, 
this time is delayed due to the small value of the mass to 
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FIG. 7: Left panel: evolution of the density perturbations Spk/p, Cfe an< A $fe towards the end of inflation with a potential given 
by Eq. (|48)l and during the ensuing period of of metric preheating. The dimensionless parameter are A ~ 1(T 12 and T = 0.05. 
The red solid line represents the Hubble radius, the red dashed line the scale i c while the blue dotted line represents the 
physical wavelength of the Fourier mode. The Fourier mode is characterized by k ~ 1.738. Right panel: same as left panel but 
with T = 0.003 and k = 1.737. The value of k is the same for the two panels but the value of k sightly differs because the end 
of inflation does not occur exactly at the same time in the two cases. 



(small in comparison with the COBE normalized value) 
and lies outside of the plot. Then, the density contrast 
grows proportional to the scale factor as discussed in Sec- 
tion nni 

V. CONCLUSION 

In this paper we have investigated the evolution of sub- 
horizon perturbations between the end of cosmic inflation 
and the beginning of a radiation dominated early Uni- 
verse. This study assumes that energy density of the 
Universe during this phase is dominated by the oscilla- 
tions of a single scalar field around its minimum in a 
quadratic or quartic potential. 

We have found that perturbations are subject to a 
preheating instability in the narrow resonance regime, 
leading to the growth of the density contrast Spk/p on 
sub-horizon scales. In the case of a quadratic poten- 
tial, the domain of instability (in terms of comoving 
wavenumbers) widens with time. This implies that a 
wavenumber that has entered the band at some time 
remains inside for the whole duration of the preheat- 
ing era. The density contrast of sub-horizon unstable 
modes then increases as the scale factor, reaching non- 
linearity for a wide range of scales. Maximum growth 
occurs for spatial scales of the order of the Hubble radius 
at the end of inflation, to which is associated a mass scale 
M ond ~ 3.3 x 10 3 g (iWl0 13 GeVr\ 



A quartic potential leads to a different phenomenol- 
ogy, although the main effect of preheating instability 
remains. Notably, the range of unstable wavenumbers is 
now fixed in terms of comoving wavenumbers and peaked 
around the Hubble scale at the end of inflation. On the 
other hand, the density contrast of sub-horizon unstable 
modes increases exponentially fast. 

The growth of structure on sub-horizon scales is known 
to occur in a matter dominated era, hence the growth 
of the density contrast for a quadratic inflaton potential 
can thus be well understood in this respect. However, 
the result obtained for a quartic potential confirms that 
the effect goes beyond a simple analogy with Newtonian 
collapse of a self-gravitating dust cloud, as the effective 
equation of state in this latter case is that of radiation. 

These results open interesting avenues of research in 
the search for observational signatures of the dynamics 
of the very early Universe. In particular, as we have 
noted, the abundance of black holes of mass comparable 
to the above is constrained by cosmological arguments. 
The formation of inflaton halos may also lead to copi- 
ous production of gravitational waves, as discussed in a 
companion paper. Finally, the very existence of non- 
linearities on small scales casts some shadow on the reli- 
ability of the predictions for curvature perturbations on 
those large scales responsible for the present large-scale 
structure. These latter could be modified if effects of 
mode-mode coupling become important. 
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